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Is the core of M87 the source of its TeV emission? Imphcations 

for unified schemes. 

Markos Georganopoulos^'^ Eric S. Perlman ^, Demosthenes Kazanas^ 

ABSTRACT 

M87 has been recently shown to be a TeV source which is likely to be variable. 
Based on this, and on contemporaneous optical and X-ray monitoring, we argue 
that the source of the TeV emission is the core of M87 and not one of two jet 
knots (HST-1 and A) with X-ray brightness comparable to that of the core. 
We model the TeV emission in the core as inverse Compton (IC) emission from 
the base of the jet. Homogeneous models fail to reproduce the spectral energy 
distribution (SED) and, in particular, the TeV flux. They also fail to comply 
with the unified scheme of BL Lacertae (BL) objects and FR I radio galaxies. A 
jet that decelerates from a Lorentz factor F ~ 20 down to F ~ 5 over a length 
of ~ 0.1 pc reproduces the observed SED of the M87 core, and, when aligned to 
the line of sight, produces a SED similar to those of TeV BLs. The TeV flux in 
the decelerating jet model is successfully reproduced as upstream Compton (UC) 
scattering, a recently identified emission mechanism, in which energetic electrons 
of the upstream faster flow upscatter the low energy photons produced in the 
slower downstream part of the flow. 

Subject headings: galaxies: active — quasars: general — quasars: individual 
(M87) — radiation mechanisms: nonthermal — X-rays: galaxies 



1. Introduction 

Gamma-ray emission in GeV energies has been detected in several bright BLs (Hartman 
et al. 1999). A handful of nearby, bright BLs have been also detected in the TeV regime (e.g. 
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Aharonian et al. 2005). Gamma-ray emission at some level is expected from radio galaxies, 
since, according to the unified scheme (e.g. Urry & Padovani 1995), BLs are FR I radio 
galaxies (Fanaroff & Riley 1974) with their relativistic jets oriented toward the observer. 
So far, the only radio galaxies detected in Gamma-rays are the FR I's Cen A (Sreekumar et 
al. 1999) and M87 {HEGRA, Aharonian et al. 2003; HESS, Beilicke et al. 2005) at TcV 
energies. Also, a tentative GeV detection of NGC 6251 (Mukherjee et al. 2002) has been 
strengthened by its hard X-ray detection. (Foschini et al 2005). 

M87 is among the nearest galaxies with a bright radio/optical/X-ray jet. Due to its 
proximity (distance 16 Mpc) the jet has been studied with unparalleled spatial resolution 
(e.g., Perlman et al. 1999, 2001 and references therein). The synchrotron emission extends 
to X-ray energies in all jet knots, as well as the core (e.g., Perlman et al. 2001, Marshall 
et al. 2002, Wilson & Yang 2002, Perlman & Wilson 2005). HST observations of the core 
(Tsvetanov et al. 1998) showed variability on timescales of ~ 1 month. HST and Chandra 
semi-monthly monitoring of the jet during the last three years (Harris et al. 2003, 2005; 
Perlman et al. 2003) established the broadband nature of this variability. Surprisingly, these 
observations detected an enormous (factor ~ 50) ongoing flare in the jet component HST-1, 
0.8" (60 pc projected) from the core, with similar variability timescales. 

The discovery of flaring optical and X-ray emission established an important link be- 
tween M87 and the BLs, as did the detection of apparent superluminal motion in the jet 
at radio (Biretta, Zhou, & Owen 1995) and optical (Biretta, Sparks, & Macchctto 1999) 
frequencies. The apparent superluminal velocities in the optical of Uapp ~ 6 c require a jet 
orientation angle 6 < 19° and F > 6. The detection of TeV emission constitutes another 
link between FR I's and BLs. The TeV data appear to indicate significant variabihty, as the 
1998/99 (HEGRA) and 2003/04 (HESS) fluxes differ by over a factor 3 at > Scr significance 
(Beihcke et al. 2005). In BLs, correlated X-ray and TeV variations (e.g. Takahashi et al. 
1996) suggest that the TeV emission is produced at the core. In the case of M87, it is nec- 
essary to consider as possible TeV emitters three components with comparable optical and 
X-ray synchrotron luminosity: the core, HST-1, and knot A, an extended jet feture located 
12.4" from the core. 

Based on the likely TeV variability, we argue in §2 that the most plausible source of 
the TeV emission is the core and not HST-1 or knot A. In §3 we consider the modeling 
constraints imposed by the core SED and the unified scheme and show that homogeneous 
models cannot satisfy these constraints. We then show (§3.1) that a decelerating relativistic 
jet can solve these problems. Finally, in §4 we sum up our conclusions. 
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2. The source of the TeV emission: Variabihty constraints. 

The production of TeV photons requires TeV energy electrons, which also emit optical - 
X-ray synchrotron radiation, resulting in correlated TeV and optical - X-ray variability. The 
possible TeV emission sites are, therefore, the brightest optical - X-ray features of the M87 
jet. These are the core and knots A and HST-1, all of comparable optical and X-ray fluxes 
(within a factor ~ 5 — 10). We now use the fact that the TeV emission has hkely dropped 
by a factor of ~ 3 between 1998/99 and 2003/04, to constrain the actual TeV source. 

We first consider knot A, one of the largest bright regions in the M87 jet, with its fiux 
maximum region Ri 1" x 0.7" in HST (Sparks, Biretta, & Macchetto 1996) and Chandra 
(Perlman & Wilson 2005) images, translating to 80 x 55 pc (projected). Such a large region 
cannot vary on timescales of a few years, and, indeed, monitoring by both ROSAT (Harris, 
Biretta, & Junor 1999) and Chandra (Harris et al. 2005) does not detect variabihty of 
more than ~ 10 — 20% in the X-rays. 

We turn now to knot HST-1, which is currently exceeding the core fiux in the optical, 
UV, and X-rays by factors > 2 — 5, depending on the band (Harris et al. 2003, 2005; 
Perlman et al. 2003). In Spring 2004, the X-ray fiux of HST-1 was ~ 3 times brighter than 
it was in Spring 2003, and since 1998/99 the increase is far larger - nearly a factor 50. By 
comparison, there is no evidence of an increase in TeV flux from 2003 to 2004, or from 1998 
to 2004 (Aharonian et al. 2003; Beilicke et al. 2005). Since the X-ray and TeV variations 
are not correlated, the contribution of HST-1 to the TeV output of M87 seems not to be 
significant. This may indicate a low synchrotron photon compactness for HST-1, which can 
be used to set, in a manner similar to that of Stawartz et al. (2004) for knot A, a lower limit 
on its magnetic field. 

Having excluded knots A and HST-1, we turn to the core. HST, ROSAT and Chandra 
monitoring (Harris et al. 1999, 2003, priv. comm.; Perlman et al. 2003) suggest that its 
variations may be similar to those observed by HEGRA and HESS in 1998/99, and 2003/04. 
We therefore believe that the most hkely source for the TeV emission is the core of M87. 

3. Modeling the core SED 

Given the ~ 1 month variability timescale of the core emission, its modeling requires 
contemporaneous broadband observations. The most complete coverage of the core SED is 
that of April and May 2003 (Figure 1). We include non-contemporaneous radio and IR data 
(crosses), given that the amphtude of the core variabihty does not exceed a factor of ~ 2 
amphtude. We also plot the 2a EGRET upper limit (Sreekumar et al. 1996), with the gray 
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band reaching up to a flux 4 times larger than the detection limit, as expected when the IC 
seed photons are the source's synchrotron photons: a change in the electron injection results 
respectively in linear and quadratic changes of the synchrotron and IC fluxes. 

The synchrotron peak between 10^^ — 10^^ Hz corresponds to a coohng break rather 
that to a cutoff of the electron energy distribution (EED), because this component exhibits 
a power- law shape that extends to X-ray energies (e.g. Marshall et al. 2002). This can be 
seen in the 2003 spectral data (Perlman et al. 2003; Harris et al. 2003). The decrease in 
the synchrotron SED from April to May 2003 must be accompanied by a quadratic decrease 
in the IC SED. If there was no synchrotron contribution in the X-rays, the X-ray spectral 
slope should have remained constant, with only a reduction of the X-ray flux. However, 
the X-ray slope steepens, suggesting a significant synchrotron contribution in the low state. 
Note that an extrapolation of the two higher energy fiuxes (1.4 and 4.0 KeV) down to 0.5 
KeV overproduces the 0.5 KeV fiux in both April and May 2003. This may be due to an 
intrinsic absorption component (Perlman & Wilson 2005). 

The SED of the core of M87, together with the FR I - BL unification (Urry & Padovani 
1995) argue against a homogeneous model for the core. The synchrotron SED of TeV BLs is 
characterized by a spectral break at ~ 10^^~^^ Hz, with the SED still increasing above this 

break, to reach its peak at ~ 10^^ — 10^^ Hz, beyond which it exhibits a cutoff (e.g. the 1997 
flare of MKN 501, Krawczynski et al. 2000). This is in strong contrast with the declining 
IR - optical - X-ray SED of the core of M87. Because relativistic motion preserves the SED 
shape, simply shifting it along the frequency and luminosity axes, a simple boosting cannot 
shift the SED of the core of M87 into the mold of the SED of the TeV BLs. 

Even if M87 is a special case that, if aligned, does not resemble the TeV BL Lacs, 
we still need to reproduce its SED. This turns out to be practically impossible for the 
homogeneous model. To sketch why this is the case, we will make use of four quantities 
derived or constrained by observations: the synchrotron peak frequency i^g = 10^^~^^ Hz, 
the synchrotron peak luminosity Lp ^ 10^^ erg s~^ , the Compton dominance (ratio of IC 
to synchrotron luminosity) ~ 1 — few, and the variability timescale t^ar ~ 1 month. The 
bolometric synchrotron luminosity is Ls^boi — Lpkb, where kb — 2/[(3 — p)(p — 2)] and p is the 
index of the injected power law electron distribution. If 7^ is the Lorentz factor where the 
cooling time equals the escape time, 75 oc B^^R~^, where B is the magnetic field and R is the 
radius of the source, and Us oc SB'-/'^ oc SB^^R^"^, where S is the usual Doppler factor. Noting 
that Cd = Us/Ub, where Ug oc Lg^^oiR '^^ '^ is the synchrotron photon energy density and 
Ub = -B^/(87r) is the magnetic field energy density, we have oc Lpk^R^^ B^'^5^^ . Using 
these expressions for Vg and C^, together with R ~ ctyarS, we obtain B oc d^Cd^J^ L^^k^^ and 
5^ oc VsLy^ k^^"^ C^^^'^t~^^. The last expression permits only a very narrow range 5 ~ 1 — 2, 
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even after allowing for a generous range in the SED-derived parameters. Using the above and 
noting that the peak of the synchrotron self Compton (SSC) component is at a frequency 
^ssc ~ ^sll, we obtain 



where 1/5,14 = z/s/lO-*^^, Lp^n — L^/IO^^, and tvar,m is the variability timescale in months. For 
a SED - derived set of parameters (I'p,4i — 1, — 1, k}, — 8.33 for p — 2.6 required by 
spectral fitting), the SSC component cannot peak at vssc ^ 10^^ Hz, even for Vg — 10^^ Hz. 

While at vssc the emission is in the Thomson regime, due to the reduced Klein- 
Nishina cross-section the SSC luminosity Lssc,kn at higher frequencies drops much faster, 
Lssc,KN oc (Georganopoulos et al., in prep), than anticipated in the Thomson case, 

Lssc,T oc i/~(p~2)/2 Because of this steeper SED at the upper end of the SSC component 
and the upper limit on pssc (see Eq. (1) above), homogeneous models fail to match the 
observed TeV flux. A typical example of this is shown in Figure 1, where we plot the SED 
for two different jet powers. Note how in the low state the stronger reduction of the SSC 
SED results in a steeper X-ray spectrum, in agreement with the observations. 



It has been noted (Georganopoulos & Kazanas 2003b, GK03b) that homogeneous models 
can fit the SED of TeV BLs only by invoking very high Doppler factors {6 ~ 50; e.g. 
Krawczynski, Coppi, & Aharonian 2002). Similarly fast fiows (Perlman et al. 2003) are 
required by homogeneous models if one wants to relate the M87 variability timescale (~ 1 
month) with those of TeV blazars (~ few hours) . These values are in strong disagreement 
with the BL - FR I unification, because they severely underpredict the observed FR I 
luminosities and overpredict their number counts (Chiaberge et al. 2000; Trussoni et al. 
2003). They also confiict with the slow apparent motions {uapp < c) in the pc-scale jets of 
TeV BLs (Marscher 1999; Finer & Edwards 2004, 2005). This motivated GK03b to propose 
that TeV BL jets decelerate at sub-pc scales from an initial Lorentz factor Fq ~ 10 — 20 down 
to F ~ few. In such flows, electrons are more energetic in the faster base of the flow and are 
radiatively cooled to lower energies as they advect downstream where the jet decelerates. 
Beaming, hence, is frequency dependent. In addition to this, upstream Compton (UC) 
scattering, a newly identified mechanism (Georganopoulos & Kazanas 2003a) in which the 
energetic electrons of the upstream fast fiow 'see' the seed photon field produced by the lower 
energy electrons of the slower downstream fiow boosted, increases the level and beaming of 
high energy Compton emission. 




Hz, 



(1) 



3.1. Decelerating flow and Upstream Compton scattering 
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We apply this model to the core of M87, assuming a deceleration from Fq = 20 to F = 5 
over a distance of 3 x 10^^ cm. We plot the SED as seen at an angle d — 13° for two different 
jet luminosities in Figure 2. In both cases most of the observed power comes from the slower 
part of the flow (dotted lines), although the TcV flux is UC emission from the fast base of the 
flow (broken lines). Note that, as in homogeneous models, a decrease of the jet luminosity 
(lower panel) results in a steeper - when lower X-ray SED, in agreement with observations. 

To check if, in addition to reproducing the SED, this model can reproduce the qualitative 
characteristics of TeV BLs, we plot in Figure 3 in thick lines the same SED viewed at 
9 = 1/Fo = 2.9° (for comparison wc rc-plot in thin lines the SED viewed at = 13°). 
While at 6' = 13°, Vs ~ lO^^'^ Hz corresponds to the EED cooling break, at 6^ = 2.9° the 
increased beaming in the fast base of the flow dominates the high energy emission, resulting 
in Vg ~ 10^^ Hz. The IC flux increases more than the synchrotron, due mostly to UC emission 
from the fast base of the flow, resulting in a Compton dominance Ca ~ 20, comparable to 
that of lES 1426+428, the most Compton dominated TeV BL (Costamante et al. 2003). 

4. Discussion 

The 1998/99 TeV flux of M87 is higher than the 2003/04 flux by over a factor 3 at 
> 3(7 signiflcance (Beihcke et al. 2005). Based on this and on optical and X-ray variability 
constraints, we argued that the source of the TeV emission of M87 is its core, as has been 
previously predicted both in the context of leptonic (Bai & Lee 2001) and hadronic models 
(Protheroc, Donea, & Reimer 2003). We excluded knot A because of its large size and weak 
or absent variability, and HST 1 because its optical and X-ray flux increased, while the TeV 
(HESS) flux remained constant or slightly decreased between 2003 and 2004. TeV emission 
from either the extended jet of M87 (Stawarz et al. 2004) or from hadronic processes in the 
giant elliptical galaxy itself (Pfrommer & EnBhn 2003) also do not agree with the observed 
TeV variabihty. Stawarz et al. (2005) used the fact that the TeV emission varies to set an 
upper hmit on the actual TeV emission of knot A (lower that the lowest TeV flux observed) , 
and, through this, a lower limit on its magnetic fleld. 

A homogeneous core model is not favored, because it cannot reproduce the core SED and 
is inconsistent with the unifled scheme. A decelerating relativistic flow (GK03b) satisfles both 
the SED modehng and uniflcation requirements. This is in contrast to both Cen A and NGC 
6251, that have be modeled as homogeneous sources (Chiaberge, Capetti, & Celotti 2001; 
Chiaberge et al. 2003). Neither of these sources is superluminal, suggesting a jet inclination 
larger than that of M87 and, in the context of a decelerating flow, a SED dominated by the 
slower part of the jet. Since this is characterized by small velocity gradients, their SEDs can 
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be accommodated by homogeneous models. An orientation difference also explains the lack 
of TeV emission from these sources: their TeV emission is more tightly beamed along the 
jet axis and points away from our line of sight. 

The postulated deceleration could be due to entrainment of external matter (e.g. Bow- 
man, Leahy, & Komissarov 1996) or Compton drag (e.g. Meha & Konigl 1989) on a photon 
field external to the high energy emitting jet, such as in the fast spine - slow sheath model 
of Ghisellini, Tavecchio, & Chiabcrge (2005). Both longitudinal (GK03b) and lateral (Ghis- 
ellini et al. 2005) velocity gradients can be present and quantifying their importance will 
require detailed variability modeling. We note that arguments for further deceleration have 
been presented for FR I jets at ~ 1 — 10 Kpc scales (Laing & Bridle 2002), and for powerful 
FR II and quasars jets at ~ lOO's of Kpc (Georganopoulos & Kazanas 2004). 

We end with a cautionary note: if M87 were at a distance similar to that of TeV BLs, 
the core and knot HST-1 would be seen as a single unresolved source in both optical and 
X-ray energies. Modeling the broadband SED as being produced at a single site could then 
result to unrealistic jet descriptions. 
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Fig. 1. — Open symbols: the core SED in April 2003; optical (Perlman et al. 2003), X- 
ray (Harris et al. 2003; Perlman et al. 2003) and TeV (Beilicke et al. 2005). Solid 
symbols: early May 2003 SED; (Perlman et al. 2003; Harris et al. 2003). Crosses: non- 
contemporaneous radio (Biretta, Stern, & Harris 1991) and 10 nm (Perlman et al. 2001) 
data (the gray bands indicating flux variations by a factor of 2). We also plot the 2a EGRET 
upper limit (Sreekumar et al. 1996), with the gray band reaching up to a flux 4 times larger 
than the detection hmit. The dotted (synchrotron and SSC) and sohd (total) hnes correspond 
to the model SEDs of a homogeneous model resulting for R = lO^^cm, B = 0.03 G, 5 = 1.4, 
derived from the SED using the formalism of §3. The two SEDs differing only in the jet 
comoving power Lj^t, with Lj^t = 7.5 x 10^^ erg for the lower one and Lj^t = 10^"^ erg 
for the higher one. The broken line indicates the analytical slope of the SSC spectrum 
in the KN regime. 
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Fig. 2. — Upper panel: The SED of a decelerating flow seen under an angle 9 = 13°. The 
flow decelerated from Fq = 20 to F = 5 in a distance 2; = 3 x 10^^ cm, following a F oc z""^ 
profile. The inlet has a radius of 5 x lO^^cm. The injected EED is a power law with slope 
p — 2, the magnetic field at the inlet isS = 0.015G, and the jet power is Ljet = 2.2 x 10^^ 
erg s^^. The solid lines represent the synchrotron, Compton, and total luminosity. The 
broken lines represent the synchrotron and Compton luminosity of the fastest fifth of the 
flow, while the dotted lines that of the rest of it. Lower panel: The same with reduced jet 
power, Ljet = 1-6 x 10"^"^ erg s~^. In both panels the data are the same as in Figure 1. 
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Fig. 3. — The SED of the same decelerating flow as in the upper panel of Figure 2, as seen 
under an angle of ^ = 13° (thin lines) and 6 = I/Fq = 2.9° (thick lines). As in Figure 2, the 
broken lines represent the synchrotron and Compton luminosity of the fastest fifth of the 
flow, while the dotted hnes that of the rest of it. Note how the fast part of the flow dominates 
the total luminosity of both the synchrotron and Compton components at ^ = I/Fq = 2.9°. 
The data are the same as in Figure 1. 



